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Abstract 

In 2004, the International Agency for Research on Cancer (lARC) reclassified formaldehyde (FA) from a probable 
(Group 2A) to a known human carcinogen (Group 1) citing results for nasopharyngeal cancer (NPC) mortality from 
the follow-up through 1994 of the National Cancer Institute formaldehyde cohort study. To the contrary, in 2012, 
the Committee for Risk Assessment of the European Chemicals Agency disagreed with the proposal to classify FA 
as a known human carcinogen (Care. lA), proposing a lower but still protective category, namely as a substance 
which is presumed to have carcinogenic potential for humans (Care. IB). Thus, U.S. and European regulatory 
agencies currently disagree about the potential human carcinogenicity of FA. 

In 2013, the National Cancer Institute reported results from their follow-up through 2004 of the formaldehyde 
cohort and concluded that the results continue to suggest a link between FA exposure and NPC. We discuss in 
this commentary why we believe that this interpretation is neither consistent with the available data from the most 
recent update of the National Cancer Institute cohort study nor with other research findings from that cohort, other 
large cohort studies and the series of publications by some of the current authors, including an independent study 
of one of the National Cancer Institute's study plants. 

Another serious concern relates to the incorrectness of the data from the follow-up through 1994 of the National 
Cancer Institute study stemming from incomplete mortality ascertainment. While these data were corrected by 
the National Cancer Institute in subsequent supplemental publications, incorrect data from the original publications 
have been cited extensively in recent causal evaluations of FA, including lARC. We conclude that the NCI 
publications that contain incorrect data from the incomplete 1994 mortality follow-up should be retracted entirely 
or corrected via published errata in the corresponding journals, and efforts should be made to re-analyze data from 
the 2004 follow-up of the NCI cohort study. 

Keywords: Formaldehyde, Nasopharyngeal cancer. Cohort mortality study. National Cancer Institute, Reanalyses, 
Silversmithing 



Introduction 

In 2004, the International Agency for Research on Cancer 
(lARC) reclassified formaldehyde (FA) from a probable 
(Group 2A) [1] to a known human carcinogen (Group 1) 
[2] citing the National Cancer Institute (NCI) cohort ana- 
lyses for nasopharyngeal cancer mortality (NPC) [3]. Based 
on the same NCI findings, the Group 1 classification was 
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upheld by lARC following the working group meeting for 
lARC Monograph Volume 100 F [4]. In Europe, the NCI 
FA cohort study as published by Hauptmann et al. [3] and 
evaluated by lARC 2006 [2] was also debated in detail to 
justify an appropriate cancer classification for FA according 
to the EU rules. The Committee for Risk Assessment of 
the European Chemicals Agency (RAC) disagreed with the 
proposal to classify FA as a known human carcinogen 
(Care. lA), proposing a lower but still protective category, 
namely as a substance which is presumed to have carcino- 
genic potential for humans (Care. IB). In 2011, in the US, 
the National Toxicology Program classified FA as a "known 
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to be a human carcinogen" based in large part on the find- 
ings of the NCI study [5]. 

The French Agency for Food, Environmental and Oc- 
cupational Health & Safety (ANSES)^ on behalf of the 
French Ministries of Health, the Environment and Labor, 
prepared a proposal for harmonized classification and 
labeling of FA*' after lARC's decision to classify as a 
human carcinogen. The responsible Committee for Risk 
Assessment (RAC)'^ of the European Chemicals Agency'* 
evaluated this proposal in 2012. The major discussion 
issue at RAC was the correct interpretation of the NCI 
FA cohort study on NPC mortality [3] taking into ac- 
count all peer-reviewed publications that referred to this 
study and that explored the issue of NPC risk after FA 
exposure, in particular those publications that were pub- 
lished after 2004. 

The Committee for Risk Assessment disagreed with 
the proposal to classify FA as known to have carcino- 
genic potential to humans (Care. lA), proposing a lower 
but still protective category, namely as a substance 
which is presumed to have carcinogenic potential for 
humans (Care. IB)*^. Thus, U.S. and European regulatory 
agencies currently disagree about the potential human 
carcinogenicity of FA. Despite the fact that the classifica- 
tion systems of the different agencies are not identical 
the main message is clear: The RAC did not follow 
lARC's evaluation on human carcinogenicity. 

The National Cancer Institute formaldehyde cohort study 

In June 2013, the NCI published the findings of its update 
through 2004 of mortality from solid tumors among 
workers in the US industry-wide FA study [6]. This study 
includes 10 plants and represents the largest cohort study 
of workers with potential exposure to FA [7]. The purpose 
of the Beane Freeman et al. update [6] was to extend the 
mortality follow-up through 2004 and to examine the as- 
sociations among different exposure characterizations and 
mortality from several solid tumors. This study also in- 
cluded corrections to the earlier update of mortality 
through 1994 published in 2004 by Hauptmann et al. 
[8]. Beane Freeman et al. [6] claim that a persistent in- 
creased risk remains for NPC mortality associated with 
peak, average intensity and cumulative FA exposure 
metrics as reported in Hauptmann et al. [3], although 
this NPC risk was not reported by Blair et al. [7] in the 
original FA cohort analysis based on follow-up through 
1979. The main conclusion from Beane Freeman et al. 
[6] is that the update through 2004 suggests a link be- 
tween FA exposure and NPC mortality that is consistent 
with some case-control studies [9-14]. Aside from not 
statistically significantly elevated rate ratios for salivary 
gland cancer mortality, the authors observed no associa- 
tions with mortality from other cancer types reported in 



other studies, including lung, laryngeal, nasal sinus and 
brain [1,4]. 

The NCI FA cohort study has a number of strengths 
including size, long follow-up time, several individual 
worker-level FA exposure metrics, partial adjustment for 
potential confounding factors, and the use of external 
and internal comparison populations. However, we have 
several concerns regarding the interpretation of the epi- 
demiological findings from this study by Beane Freeman 
et al. [6], in particular, the purported link between FA 
exposure and NPC mortality. We believe that this inter- 
pretation is neither consistent with the available data 
nor with other research findings based on this group of 
US formaldehyde workers. The series of publications by 
the current authors and others present evidence [15-19]. 
Another serious concern relates to the incorrectness of 
the data reported by Hauptmann et al. [3,20] stemming 
from incomplete mortality follow-up that was first re- 
ported by Beane Freeman et al. [21]. While these data 
were corrected in subsequent supplemental publications 
[8,22], incorrect data from the original publications by 
Hauptmann et al. have been cited extensively in recent 
causal evaluations of FA, including lARC [4]. We note 
that lARC's evaluation in 2004 [2] was mainly based on 
Hauptmann et al. [3,20]. 

This commentary describes our concerns and provides 
justification for continued re-analysis of the NCI cohort 
data on mortality. We will refer to three specific publica- 
tions from updates of the NCI FA cohort: namely, the 
1994 mortality update [3], corrections to the 1994 mor- 
tality update published as supplementary data [8,22] and 
the 2004 mortality update [6]. 

Issue 1 - Inconsistently reported cohort data 

In 2010, Marsh et al. [15] described several major discrep- 
ancies between the number of observed deaths reported 
by Hauptmann et al. [3] and Beane Freeman et al. [22]. A 
main finding was that the corrected increase for total 
deaths by FA exposure status apparentiy equals 995, not 
1006 as described by Beane Freeman [21]. We also found 
that when exposure is defined using a 2-year lag (used in 
the NCI analysis of lympho-hematopoietic cancers [3]), 
the percent increase in corrected numbers of deaths 
among FA "unexposed" subjects is approximately two 
times greater than that observed among FA "exposed" 
subjects for all deaths, all cancer deaths and all solid can- 
cers [15]. Moreover, the basis of NCI's revised count of 
1006 deaths is neither evident nor reproducible from in- 
formation provided. The most recent NCI publication 
by Beane Freeman et al. [6,8] did not clarify the discrep- 
ancy. Beane Freeman et al. [6] state (page 1016) that 
Hauptmann et al. [3] had missed 999 deaths and had 
counted four living subjects as deceased. They also re- 
ported 20 subjects whose age had exceeded the age of 



Marsh et al. Journal of Occupational Medicine and Toxicology 2014, 9:22 
http://www.occup-med.eom/content/9/1/22 



Page 3 of 9 



the oldest known decedent (102 years) as deceased. This 
information suggests that the corrected total death 
count should be 1015 (999-4 + 20), not the 1006 re- 
ported by Beane Freeman et al. [6]. Inconsistent with 
both numbers (1015 and 1006) that can be derived from 
Beane Freeman et al. [6], Table SI in Beane Freeman 
et al. [8] yields a difference of 995 in comparison to 
Table two in Hauptmann et al. [3] . The identified differ- 
ences are small and may have no substantial effect on 
estimates of relative risks, but the discrepancies between 
the original and corrected total numbers of deaths in 
the NCI cohort creates doubt whether all ascertainment 
errors were addressed by the NCI, as the consistent 
differential change in cause-specific numbers of deaths 
by exposure category cannot likely be explained as a 
chance occurrence. At the time of this publication, the 
NCI investigators have not presented an explicit clarifi- 
cation of the mortality ascertainment errors and the 
changes in total deaths and cause-specific deaths from 
the reported numbers presented in the Hauptmann 
et al. publications. 

Issue 2 - Relevance of NCI's corrected number of deaths 

Tables 1, 2, 3 summarize key results from Hauptmann 
et al. [3], the only published corrections to this publica- 
tion [8], and the updated mortality results reported by 
Beane Freeman et al. [6] for all causes of death com- 
bined and for NPC mortality in relation to FA exposure 
status based on external and internal comparisons. The 
inclusion of the claimed 1006 deaths missing in the ana- 
lyses by Hauptmann et al. [3] had a major impact on the 
corrected risk estimates for total mortality, especially 
those based on internal cohort comparisons (Table 2). 
For example, excesses in total mortality ranging from 7% 



to 21% among workers in the two highest categories of 
ever-peak FA exposure (2-4 ppm and > 4 ppm) attenu- 
ated to deficits ranging from 4% to 15% in analyses that 
corrected the Hauptmann et al. publication [8] or updated 
the mortality follow-up [6], including some reduced risk 
estimates that were statistically significant. The overall 
standard mortality ratios (SMRs) were relatively unaffected 
by the corrections but corrected results for the 1994 mor- 
tality follow-up were more consistent with those reported 
in the most recent 2004 follow-up (Table 1). Thus, to clar- 
ify the record on the impact of these corrections, we sub- 
mit that the incorrect NCI publications [3,20] should be 
retracted entirely or corrected via published errata in the 
corresponding journals. 

Issue 3 - New NCI findings for NPC argue against an 
association with FA exposure 

Beane Freeman et al. [6] reported results for updated 
mortality ascertainment through 2004 that provide evi- 
dence against an association between FA and NPC mor- 
tality. First, the most recent 10-year period of follow-up 
of the NCI FA cohort yielded only one additional NPC 
death. Further, this single NPC death occurred in the 
lowest exposure categories for the three FA exposure 
metrics evaluated by Beane Freeman et al. (Table 3). The 
observation of only one NPC death during a 10-year up- 
date of this large cohort is less than the corresponding 
number of expected NPC deaths for this period (1.14 as 
derived from information provided in the supplemental 
tables and the published report by Beane Freeman et al. 
[6,8]). First, this result by itself, argues against an associ- 
ation between FA exposure and NPC mortality. Second, 
unlike results reported in Hauptmann et al. [3] and the 
corrected analyses for the update through 1994 [8], the 



Table 1 Standardized mortality ratios and numbers of deaths for all causes combined by NCI update and FA exposure 
status^ 



NCI update 




Formaldehyde exposure status 








Non-exposed 




Exposed 




No. deaths 


SMR 95% CI 


No. deaths 


SMR 


95% CI 






All causes of death (ICD8: 001-999) 








NCI 1 994 update"- 


1,991 


0.85* 0.81-0.89 


6,495 


0.96* 


0.94-0.98 


NCI 1 994 update (corrected)''- 


2,169 


0.89* 0.86-0.93 


7,312 


1.02 


0.99-1.04 


NCI 2004 update"- 


2,605 


0.90* 0.87-0.94 


11,346 


1.03* 


1.01-1.05 






NPC deaths (ICD8: 147)^^- 








NCI 1 994 update''- 


2 


1 .56 0.39-6.23 


8 


2.10 


0.91-414 


NCI 1 994 update (corrected)''- 


2 


1 .57 040-6.28 


8 


2.13 


0.92-419 


NCI 2004 update=- 


2 


1.45 0.17-5.25 


9 


1.84 


0.84-3.49 



^ Hauptmann et al. [3]. 
"^ Beane-Freeman et al. [8]. 
"^ Beane-Freeman et al. [6]. 

■^ 95% confidence intervals based on exact methods. 
^"FA exposure status calculated using 15-year time lag. 
*95% CI does not include 1.00. 
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Table 2 Rate ratios and numbers of deaths for all causes combined by NCI update and peak FA exposure' 



Peak exposure (ppm) 







0 




>0 - < 2.0 


2.0 - < 4.0 


>4.0 








NCI update 


RR 95% 


CI 


No. 
deaths 


RR 95% CI 


No. 
deaths 


RR 95% CI 


No. 
deaths 


RR 95% CI 


No. 
deaths 


p-value 
trend'' 


p-value 
trend'- 


NCI 1 994 update"- 


1.05 




1,991 


1.00 Baseline 


2,554 


1.21* 


1,945 


1.07* 


1,996 


0.013 


0.014 


NCI 1 994 update 
(corrected)''' 


1.02 0.85- 


1.04 


2,169 


1.00 Baseline 


3,201 


0.96 0.91-1.02 


2,012 


0.86*0.81-0.91 


2,099 


<0.001 


<0.001 


NCI 2004 update^- 


0.98 0.93- 


1.04 


2,605 


1.00 Baseline 


4996 


0.95* 0.90-0.99 


3,096 


0.85*0.81-0.89 


3,254 


<0.001 


<0.001 



" Hauptmann et al. [201. 
"^ Beane-Freeman et al. [81. 
"^ Beane-Freeman et al. [61. 

^ Two-sided likelihood ratio test (Idf) of zero slope for continuous formaldehyde exposure among exposed person-years only. 
''"Two-sided likelihood ratio test (Idf) of zero slope for continuous formaldehyde exposure among unexposed and exposed person-years. 
^FA exposure metric calculated using 15-year time lag. 
*95% CI does not include 1.00. 



overall SMR for NPC among all FA-exposed workers in 
the most recent update through 2004 [6] is no longer 
greater than 2.00 and is not statistically significant 
(SMR = 1.84, 95% CI = 0.84-3.49) (Table 1). Third, the 
SMR for NPC among FA-exposed workers is now con- 
sistent with the corresponding SMR for all workers 
unexposed to FA (SMR = 1.45, 95% CI = 0.17-5.25), as 
the 95% confidence interval for FA-exposed workers is 



entirely contained within the confidence interval for 
FA-unexposed workers (Table 1). 

Table 3 also shows that for the ever-peak exposure 
metric, the rate ratio (RR) for the two NPC deaths in the 
unexposed category compared with the single NPC death 
in the lowest exposed category is now 4.39 (95% CI = 0.36- 
54.05). This result is consistent with the RR for the highest 
category of ever-peak exposure (RR = 7.66, 95% CI = 0.94- 



Table 3 Rate ratios and numbers of deaths for NPC by NCI update and FA exposure metric' 



NCI update 



RR 95% CI 



No. RR 95% CI No. 
deaths deaths 



RR 95% CI 



No. 
deaths 



RR 95% CI 



No. p-value 
deaths trend** 



p-value 
trend'- 



0 

1.00 Baseline 
1.00 Baseline 
NCI 2004 update^- 439 0.36-54.05 



NCI 1 994 update"- 

NCI 1 994 update 
(corrected)'^' 



0 

NCI 1994 update" 1.00 Baseline 

NCI 1 994 update , „ , 

, |,f 1 .00 Baseline 

(corrected) 

NCI 2004 update"- 6.79 0.55-83.64 



NCI 1 994 update"- 240 

NCI 1994 update 2,41 0.35-16.70 
(corrected) - 

NCI 2004 update"- 1 .87 0.30-1 1 .67 



Peak exposure (ppm) 
> 0 - < 2.0 2.0 - < 4.0 

NA 0 NA 0 



NA 



0 



NA 



1 .00 Baseline 1 NA 

Average Intensity of exposure (ppm) 
0.1-04 0.5-0.9 
NA 0 1.38 



NA 



0 



0.37 0.03-4.63 



1 .00 Baseline 1 2.44 0.15-39.07 

Cumulative exposure (ppm-years) 
>0-<1.5 1.5 -< 5.5 

1 .00 Baseline 3 1.19 



1.00 Baseline 
1.00 Baseline 



1.20 0.12-11.56 
0.86 0.10-7.70 



4.0 

1.83 

1 .82 0.32-1 0.46 
7.66 0.94-62.34 

> 1.0 

1.67 

1 .66 0.29-9.48 
11.54* 1.38-96.81 

> 5.5 

414 
4.15 0.83-20.78 
2.94 0.65-13.28 



0.044 

<0.001 
0.005 

0.126 
0.07 
0.09 

0.029 
0.04 
0.06 



<0.001 
0.05 
0.10 

0.066 
0.14 
0.16 

0.025 
0.05 
0.07 



^ Hauptmann et al. [31. 
"^ Beane-Freeman et al. [8]. 
'^ Beane-Freeman et al. [61. 

''Two-sided likelihood ratio test (Idf) of zero slope for continuous formaldehyde exposure among exposed person-years only. 
''Two-sided likelihood ratio test (Idf) of zero slope for continuous formaldehyde exposure among unexposed and exposed person-years. 
^ All FA exposure metrics calculated using 15-year time lag. 
•95% CI does not include 1.00. 
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62.34, 7 deaths), as the 95% confidence limits show a large 
overlap (heterogeneity chi-squared = 0.07 (d.f. = 1), p = 
0.79). Likewise, RRs for the unexposed categories of 
average intensity and cumulative exposures (compared 
with the lowest exposure category for both metrics as the 
referent) were also elevated (RR = 6.79, 95% CI = 0.55- 
83.64 and RR = 1.87, 95% CI = 0.30-11.67, respectively). 
Finally, for average intensity and cumulative exposures, 
RRs for the moderate exposure categories (compared with 
the lowest exposure categories for both metrics) are now 
less than the RRs for the unexposed category (RR = 2.44, 
95% CI = 0.15-39.07 and RR = 0.86, 95% CI = 0.10-7.70, 
respectively). These reported findings by Beane Freeman 
et al. [6] indicate that there is no meaningful difference in 
the RRs for NPC mortality between the unexposed and 
FA-exposed groups. 

Compared with the original and corrected 1994 mor- 
tality updates [3,8], the 2004 mortality update reported 
by Beane Freeman et al. [6] provides little or no evidence 
of an exposure-response relationship for FA exposure and 
NPC mortality as shown in Table 3. First, the reported 
RRs are difficult to compare to RRs reported earlier due to 
the substantial methodological change of using the lowest 
exposure category as the reference rather than the use of 
the unexposed category as the reference in the previous 
1994 update (discussed below under Issue 4). Second, for 
the same reason, the claim by Beane Freeman et al. [6] 
that a statistically significant trend (p = 0.005) for RRs for 
the ever-peak exposure analyses as based on FA-exposed 
workers only is difficult to compare to earlier findings (as 
discussed under Issue 5). Third, none of the trend tests in- 
corporating all exposure categories in Beane Freeman 
et al. [6] are statistically significant. As noted by Beane 
Freeman et al. [6] (page 1018), "RRs calculated using non- 
exposed as the referent were attenuated, which is also 
reflected in the tests for trend using all person-years. . . ". 
This is evident in the results summarized in Table 3 that 
p-values for all trend tests estimated based on all unex- 
posed and exposed person years are greater than 0.05. 
These results suggest that the mortality excesses for NPC 
among subjects unexposed to FA are simply due to chance 
or reflect the presence of some other exogenous factor 
outside of the study plants, consistent with findings re- 
ported by Marsh et al. [16] for Plant 1 (discussed below 
under Issue 7). 

Issue 4 - Inappropriateness of excluding unexposed workers 
from the evaluation of exposure-response relationships 

As noted throughout this commentary, Beane Freeman 
et al. [6] fails to address statistical evidence that non- 
exposure to FA is also associated with NPC mortality. 
The language used by Beane Freeman and colleagues in 
the published abstract states that previously observed ex- 
cesses of NPC mortality "persisted" in the 2004 mortality 



update. This term implies that updated analyses were con- 
sistent with the previous reports and resulted in similar 
risk estimates for the three exposure metrics: peak, aver- 
age intensity, and cumulative exposure as reported by 
Hauptmann et al. [3]. This interpretation is not accurate. 
In fact, the addition of one NPC death during the recent 
10 years of additional follow-up changed the analytic ap- 
proach for calculating RRs for both peak and average in- 
tensity exposures (Table 3). That is, the additional death 
was assigned to the "low" exposure group for each metric 
(0-2 ppm for ever-peak, and 0.1 - 0.4 ppm for average 
intensity exposure). In the previous mortality update 
through 1994 [3], and correction to this publication re- 
ported by Beane Freeman et al. [8], these "low" categor- 
ies had zero observed NPC deaths, thus the "unexposed" 
category was used as the referent for internal RR analyses. 
Therefore, the estimates reported by Beane Freeman et al. 
[6] do not represent a "persistent" excess risk for ever- 
peak and average intensity exposures as they are relative 
to different referent categories for each exposure metric 
and must be interpreted in the context of the specified ref- 
erent category. 

It is inappropriate to exclude unexposed workers from 
internal analyses as done by Beane Freeman et al. [6]. 
All workers are from the same factories and, as noted 
by McLaughlin et al. [17] in a response to a letter by 
Hauptmann and Ronckers [23], lagging of FA exposure 
by 15 years results in contributions to the unexposed 
category from workers who were, in fact, exposed to FA. 
In fact, most of the person-time at risk allocated to the un- 
exposed category represents years of follow-up of workers 
who were eventually exposed to FA. Thus, even speculated 
differences between unexposed and exposed workers with 
regard to unknown confounders cannot be used to justify 
the choice of the low-FA exposure category rather than the 
non-exposed category as the reference category for the cal- 
culation of relative risks [17]. 

Issue 5 - The trend tests used in the NCI 2004 updates 
produce misleading results and may be mis-specified 

Beane Freeman et al. [6] used trend tests for categorical 
variables based on the likelihood ratio for the slope of 
the corresponding continuous variable, except for peak 
exposure, which was based only on categorical scores. 
Beane Freeman and colleagues provide no slope param- 
eter estimates from the Poisson regression analysis of the 
continuous exposure values, and no indication that the 
authors evaluated the underlying assumptions and fit of 
the model. Assuming that the underlying assumption of 
linearity was met, Beane Freeman et al. could have used 
the slope parameter estimates to determine the increase in 
the RR per unit increase of FA exposure. In addition, 
Beane Freeman et al. placed the continuous variable-based 
trend test p-values in juxtaposition to the RRs derived 
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from the categorical analyses presented, producing seem- 
ingly inconsistent results. For example, as shown in Table 3 
the pattern of RRs across the four cumulative exposure 
categories are 1.87, 1.0 (baseline), 0.86 and 2.94. As shown 
in Figure 1, with or without the inclusion of the unex- 
posed, these data provide no evidence on an increasing 
linear trend in RRs relative to increasing cumulative ex- 
posure to FA, yet trend test p-values are reported as 0.06 
and 0.07 [6] without and with the unexposed category, re- 
spectively. Also, given that these p-values are 2-sided and 
a large number of person-years were accrued in the FA- 
unexposed category, it is also possible that these p-values 
indicate a nearly statistically significant inverse trend. 

It is also quite possible that the inconsistent pattern of 
findings between the continuous and categorical data 
reflects the presence of an underlying non-linear relation- 
ship between log RR and continuous exposure. For ex- 
ample, if the underlying log RR vs. continuous exposure 
relationship was quadratic, the continuous model includ- 
ing the unexposed subjects could yield a marginally statis- 
tically significant positive slope in conjunction with no 
apparent trend or a U-shaped relationship in the RRs 
across the four exposure categories. An informative 
Poisson regression modeling approach would present an 
evaluation of the linearity assumption of the continuous 
model (i.e., linearity of log RR vs. continuous exposure), 
and if linearity is met, the corresponding slope estimates 
and confidence intervals. 

In addition to the continuous model estimates for the 
continuous exposure variables, Beane Freeman et al. [6] 
should show trend tests based on the midpoints (or 



scores) of the categorical data presented in the tables. This 
categorization of the continuous variable is arbitrary, and 
the concurrent trend tests based on categorical midpoints 
(which also require underlying linearity with respect to log 
RR vs. midpoint score) directly reflect the reported and 
apparent pattern of categorical RRs presented in the ta- 
bles. This avoids possible misrepresentation or confusion 
about trends (or lack thereof) based on the underlying 
and unseen continuous data. Non-linearity in the continu- 
ous log RR vs. continuous FA exposure relationship, such 
as a quadratic relationship, would invalidate the current 
trend tests. Another problem with the trend tests reported 
by Beane Freeman et al. [6] is the instability of the RR esti- 
mates based on only one observed death that occur in 
some baseline categories of the reported trends. This is es- 
pecially problematic for the categorical analysis of NPC 
highest peak exposure (Table 3). 

Issue 6 - Failure to recognize the important interaction 
structure between plant group and FA exposure 

In regard to analyses reported by Marsh et al. [18], Beane 
Freeman et al. [6] state in the discussion section (page 
1023), "Based on two groups of plants (plant 1 vs. plants 
2-10) and a continuous version of the originally categorical 
peak metric they demonstrated no significant interaction 
between plant group and peak exposure (p = 0.09) and ad- 
justment for plant group did not markedly change the risk 
related to peak exposure". Beane Freeman et al. [6] further 
state, "These results are consistent with our own analyses 
that showed no evidence of plant heterogeneity for a broad 
group of metrics, including peak exposure. " First, the so- 



15 
14 
13 
12 
11 
10 

i 8^ 

i 7 
I 

: 6 
5 
4 
3 
2 
1 
0 



Trend p-value = 0.06 (exposed only) 
I Trend p-value = 0.07 (unexposed and exposed) 
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called "influence analysis" conducted by Beane Freeman 
and colleagues, amounts to excluding each of the 10 study 
plants one at a time to evaluate the consistency of the 
findings (a type of sensitivity analysis). This influence ana- 
lysis was new to the 2004 mortality update [6] and pre- 
sumably an attempt to address the heterogeneity in results 
for NPC by plant reported in the reanalyses of Marsh and 
Youk [19] and Marsh et al. [18]. Although Beane Freeman 
and colleagues address plant heterogeneity systematically, 
the authors neither correctly interpret the results of this 
analysis nor correctly interpret the results of the inter- 
action evaluation performed by Marsh et al. [18]. 

Regarding the results of the reported influence analysis 
in the publication, Beane Freeman et al. [6] conclude 
that there is an apparent absence of plant heterogeneity 
upon removing Plant 1. However, the RR for NPC mortal- 
ity in the highest category of ever-peak exposure de- 
creased 56% from 7.66 for the 7 deaths including Plant 1, 
(95% CI = 0.94-62.34) to 3.36 for 2 deaths excluding Plant 
1, (95% CI = 0.30-37.27). Moreover, the corresponding RR 
for the highest average intensity category of FA exposure 
decreased 65% from 11.54 for 6 deaths including Plant 1 
(95% CI = 1.38-96.81) to 4.09 for 1 death excluding Plant 
1, (95% CI = 0.25-66.0). These results for the influence 
analysis are quite similar to the RR for unexposed workers 
(RR = 4.39, 95% CI = 0.36-54.05). Further, when Plant 1 
was included in the analysis and any other plant (from 
Plant 2 through Plant 10) was excluded, the results were 
similar to the results for analyses reported for all 10 plants. 
This pattern of findings in Beane Freeman et al. [6] is en- 
tirely consistent with the results of Marsh and Youk [19] 
which showed that the conclusion of a causal association 
was driven heavily by the large, statistically significant ex- 
cess in NPC mortality risk for employees from Plant 1 for 
the highest ever-peak exposure category. 

Regarding the interpretation by Beane Freeman et al. 
[6] of the results of the Marsh et al. [18] interaction 
evaluation, it is not true as stated by Beane Freeman and 
colleagues (page 1023) that the Marsh et al. analysis does 
not demonstrate a significant interaction between plant 
group and peak exposure. Rather, Marsh et al. [18] found 
that after dropping the main effect of plant group from the 
full interaction model (Model 9) the estimation process 
yielded more stable findings (Model 8) and the interaction 
between the plant group indicator and continuous peak 
exposure was found to be significant (p = 0.03). Unfor- 
tunately, a rigorous interaction analysis, such as that 
recommended and conducted by Marsh et al. [18] was 
not performed by Beane Freeman and colleagues. 

Issue 7 - Misrepresented findings from tfie independent 
study of Plant 1 

In 2007, Marsh et al. [16] reported findings from a 2003 
update of their independent and expanded study of Plant 



1 designed to investigate alternative explanations for the 
anomalous findings for NPC in this single plant. The re- 
sults of the nested case-control study suggested that the 
large NPC mortality excess in Plant 1 may not be due to 
FA exposure, but rather reflects the influence of external 
employment in the ferrous and non-ferrous metal indus- 
tries of the local area that entailed possible exposures to 
several suspected risk factors for upper respiratory system 
cancer (e.g., sulfuric acid mists, mineral acid, metal dusts 
and heat) [16]. 

In response to these findings, Beane Freeman [6] state 
(page 1022), "Silversmithing was associated with risk of 
NPC, and Marsh et al. [16] concluded that the observed 
association between FA and NPC may actually be due to 
silversmithing. However, risks related to duration, cumu- 
lative and average intensity of exposure did not decrease 
when smoking and silversmithing or other metal work 
were added as adjustment variables, indicating that sil- 
versmithing does not confound the association between 
FA and NPC". We disagree with this interpretation be- 
cause it assumes the presence of an association between 
FA and NPC that was not evident in the Marsh et al. Plant 
1 study. In fact. Marsh et al. [16] observed considerable 
evidence that the overall NPC excess was not related to 
FA exposure in Plant 1. This additional evidence, sum- 
marized below, was not considered by Beane Freeman 
and colleagues in their discussion of the Marsh et al. 
Plant 1 study. 

First, large, statistically significant NPC excesses were 
observed in both short-term (less than one year) and 
long-term workers (SMR = 5.35 and 4.59, respectively). 
Second, the seven NPC cases were associated with very 
short periods of Plant 1 employment (4 of 7 NPC cases 
worked less than one year, 5 of 7 cases worked less 
than five years), with low average intensity of FA ex- 
posure (the range of exposures for the 7 cases was 
0.03-0.60 ppm with median of 0.14 ppm). Moreover, 
the NPC cases were concentrated among workers hired 
at Plant 1 during 1947-56, which was not the time 
period of highest FA exposure. Third, for average in- 
tensity of exposure, the metric most closely associated 
with the highest peak metric used in the NCI study, 
there was no evidence of an increasing NPC risk with 
increasing exposure (OR = 1.00, 11.41 and 2.18 for the 
categories <0.03 ppm (baseline), 0.03-0.159 ppm, and 
0.16+ ppm, respectively). Finally, Beane Freeman et al. 
[6], observed mostly null findings for respiratory cancer 
sites other than NPC, in contrast to elevations in risk 
for these sites in the Marsh et al., Plant 1 study. For ex- 
ample, with the exceptions of cancer of the tongue and 
gum, local county rate-based SMRs in the Marsh et al. 
study were elevated for all upper respiratory cancer sites, 
including the nasal sinuses that do not seem to be pene- 
trated by FA exposures [24] . 
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Issue 8 - Failure of other major FA studies to find an 
increased risk of NPC 

As noted by McLaughlin and Tarone in a commentary 
[25] and recent letter to the editor and response [26,27], 
Beane Freeman et al. [6] cite the NCI embalmer case- 
control study [28] in their discussion as evidence for a 
possible leukemia risk among FA-exposed workers, but ig- 
nore the absence of increased risk of NPC in the em- 
balmer study. The odds ratio for NPC among embalmers 
was 0.10 (95% CI = 0.01-1.20), despite the fact that em- 
balmers were reported to have the highest peak exposure 
to FA of any known occupation [28]. Also, the recent up- 
date of the NIOSH garment industry worker cohort [29] 
was not mentioned by Beane Freeman et al. [6] but only 
the older publication [30]. Interestingly, Meyers and 
colleagues still found no NPC deaths although follow- 
up of this FA exposed cohort was extended by 10 years 
to 2008 (number of expected deaths = 1.33). Finally, the 
most recent update of the British industry-wide FA 
study also found fewer NPC deaths than expected (1 death 
versus 2.0 expected) [31]. 

Conclusions 

The interpretation of an association between FA expos- 
ure and NPC mortality is neither consistent with the 
available data from the most recent mortality update of 
the FA cohort nor with other research findings from this 
cohort, other large cohort studies [29-31] and the series 
of publications by the current authors. The NCI publica- 
tions [3,20] that contain incorrect data from the incom- 
plete 1994 mortality follow-up should be retracted entirely 
or corrected via published errata in the corresponding 
journals. Efforts are underway by two of us (GM and PM) 
to re-analyze the cohort data from NCI's 2004 follow-up 
with focus on a rigorous evaluation of the relationship be- 
tween FA and NPC. 

Endnotes 

''http://www.anses.fr 

''Based on Regulation (EC) No 1272/2008 (CLP Regu- 
lation), Annex VI, Part 2 

■^http:// echa.europa.eu/ about-us/who-we-are/committee- 
for-risk-assessment 

■^ECHA, http://echa.europa.eu 

'^http://echa.europa.eu/web/guest/view-article/-/journal_ 
Content/c89bdbl3-09e9-497c-8e73-ddael3a842c8) 
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